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Abstract 

> 

' In hyperbolic reductions of the Einstein equations the evolution of gauge 

£Sj , conditions or constraint quantities is controlled by subsidiary systems. 

We point out a class of non-linearities in these systems which may have 
the potential of generating catastrophic growth of gauge resp. constraint 
violations in numerical calculations. 
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1 Introduction 

Most numerical calculations of solutions to Einstein's field equations are being 
plagued by an undesirably fast growth of constraint violations. A huge variety 
of reduced equations has been derived from Einstein's equations with the hope 
of finding versions with stable propagation properties and there have been sug- 
gested modifcations of the equations which were hoped to force back the solution 
to the constraint manifold (we refer to [H] for further dicussion and references). 
More recently, there have been performed stability analyses of the subsidiary 
systems which control the evolution of the gauge conditions or the constraint 
quantities (cf. pQ, |Hj and the references given there). These led to requirements 
on the coefficients of the equations which in the case of pQ resulted in geometric 
conditions on the foliation underlying the evolution by the main system. Nev- 
ertheless, the field appears to be wide open, most of the suggested remedies are 
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experimental, and the cause of the problems is not understood. In pQ, |S] the 
subsidiary systems have been considered as linear systems on given space-times. 
In contrast, we wish to emphasize in this note the non-linearity of the subsidiary 
system, which appears to have a potential of generating catastrophic constraint 
violations. To understand the extent to which these non-linearities may affect 
the numerical construction of space-times and to develop, if necessary, ways to 
avoid their effects, further investigations are needed. 

2 The hyperbolic reduction 

We shall use a few facts about the hyperbolic reduction procedure by which 
the geometric initial value problem for Einstein's field equations is reduced to 
a Cauchy problem for a hyperbolic system. There exist now many versions 
of it, whose general underlying structure are, however, more or less the same. 
Because it leads to concise expressions, we will use the metric coeffficients as 
basic unknowns and consider the representation of the field equations in which 
their evolution is governed by a system of wave equations. To emphasize the 
independence of our discussion of any particular coordinate system we shall 
employ the notion of a gauge source function introduced in [2]. We refer to [3], 
01 for more details on the reduction procedure. 

Let (M, g) denote a smooth 4-dimensional Lorentz space with smooth space- 
like Cauchy hypersurface S and U 3 x x — > F^(x x ) S V a smooth map of an 
open subset U into another subset V of K 4 . With functions x M and their (4- 
dimensional) differentials dx v prescribed on some open subset W of S one can 
solve near W the Cauchy problem for the semi-linear system of wave equations 



If the dx^ have been chosen linearly independent on W the solution will provide 
a smooth coordinate system x x on some neighbourhood of W in M. In terms 
of these coordinates the equations above take the form 



where the T M denote the contracted Christoffel symbols = g vr > T v ^ n of g^. 
This shows (ignoring subtleties arising in situations of low differentiability) that 
the contracted Christoffel symbols can locally be made to agree with any pre- 
scribed set of functions and that these function and the initial data determine 
the coordinates uniquely. We refer to these functions as gauge source junctions. 
Assume now that (M, g) is to be obtained by solving a Cauchy problem for 
Einstein's vacuum field equations. We shall derive the properties which will 
help us formulate this problem as a Cauchy problem for hyperbolic equations. 
The contracted Christoffel symbols are of particular interest to us because the 
Ricci tensor of g can be written in the form 



V„ V x" = 




r"(x A ) = -F^ l (x x ), 



R. 
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T p s „ + 2Ts\g sp g Mp r i/) %. 

(2.1) 
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Here the contracted Christoffel symbols (and the functions F v = g vpu consid- 
ered in the following) are being formally treated as if they defined a vector field 
(which, of course, they do not). Thus T v = g Vfl and V^I^ = d^Y v — x „ Y\. 
The discussion above suggests replacing in i|2.1[l the functions Y v by freely cho- 
sen gauge source functions F v (so that the resulting expression will depend in 
general on the coordinates x x not any longer only through the g^v). With this 
replacement the vacuum field equations take the form 

= = -\ g Xp g^.Xp+V^ F v) +r A " „ g vS g x » T p s „ + 2 T 5 A „ g s ? g XQl T u) \, 

(2-2) 

of a system of wave equations for the g^ v . We refer to (|2.2|) as the main evolution 
system or the reduced equations. For this system the Cauchy problem for g pv 
with data satisfying the constraints on a space-like hypersurface S is well posed. 
Suppose that g^ v is a solution of this problem near S. Since equation (|2.2|) is 
in fact of the form 

Rfiv = V( M Q v ), (2.3) 

where = — F^ with the calculated from g^ Vl it is not clear a priori 
whether the solution g will indeed satisfy the gauge condition T „ = F^ and thus 
the vacuum field equation i? M „ = 0. 

The Bianchi identity, which holds for any metric, implies 

n = 2V»(R^-^Rg^) =V M VQ„ + .R' 1 „Q M , (2.4) 

and thus a system of wave equations for the quantities Q v . We refer to this 
system as subsidiary system. 

The detailed analysis shows that if the Cauchy problem for the main evolution 
system is arranged so that the initial data satisfy the constraints and the gauge 
condition Q v = on S, it follows for the solution of the main evolution system 
that also dQ v and thus any 'time derivative' dt Q v transverse to S vanishes on 
S. The uniqueness property of the subsidiary system therefore implies that the 
solution to the main evolution system does indeed satisfy Q v = on the domain 
of dependence of S with respect to g. This reduces the local Cauchy problem 
for Einstein's field equations to the problem of solving equation 12. 2|) . which 
thus takes the central role in the analytic discussion. Of the subsidiary system 
only the homogeneity and the resulting uniqueness property arc needed. 

3 The non-linearity of the subsidiary equation 

The main evolution system is also central in numerical discussions. Because 
the initial data and dt Q v on the initial hypersurface S come with an error, 
the evolution properties of the subsidiary system will, however, also become 
important. 

There is no way to relate a numerical solution to the solution of the continuum 
problem one wants to approximate. To get some idea how errors infiltrate into 
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the various systems, it is useful to consider an analogy accessible to analytic 
methods. We assume the main evolution system to be satisfied by fields g^ u of 
class C 3 with an error term E pv of class C 1 so that 

Rfiu = V( M Q v ) + Efi,,. (3.1) 

Nothing will be assumed about the origin and structure of this error and the 

errors in the initial data Q M and dt Q v on S. 

Using the Bianchi identity with (|3.1|l gives the analogue 

V M V^ 1 Q u + R" v Q M = -2 V{E^ - 1 9fil/ E p P), (3.2) 

of the subsidiary equation l|2.4|l . This equation will in general not be homo- 
geneous any longer. There may be a way to avoid this problem. Assume a 
splitting of the form 

E^v = V( p e„) + f j 

with some C 2 vector field and some symmetric C 1 tensor field With 
l|3.1|l this gives 

Ru,u = V (/t Q' v) + f„u with Q'„ = Q u + e v . (3.3) 
If ffiv would vanish this would just amount to a redefinition 

F — > F' — F — p 

of the gauge source function which may be quite harmless. Moreover, we would 
get a homogenous wave equation for Q'^. If E^ is known but not of the form 
Vjju e^) , a homogenous system can be obtained by choosing the splitting suitably. 
If we require e M to solve the system of wave equations 

V M V e„ + R" „ e„ = 2 V f - J ^ £ p " 
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it follows from the splitting above that 



and the Bianchi identity implies for Q'^ the homogeneous equation 

V^Ql + R^Q'^0. (3.4) 

One could now try to analyse how errors in and d t Q v on S are propagated, 
in dependence of E^ v , by l|3.2|l or how errors in Q'^ and d t Q' u on S will be 
propagated by l|3.4l) . Assuming in this analysis E^ and whence i? Mt ,[g] 
as suitably bounded but given with no further information, standard energy 
estimates admit but cannot exclude exponential growth of Q M resp. Q'^. One 
might be able to work out conditions on R pu [g] which will allow one to keep the 
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errors under control. However, by itself this will be of limited use. Inserting 
l|3.1H into H3.2|) and, observing the splitting, into l|3.4|) . gives 



V M V Q v + Q A (V (A Q v) + E Xv ) = -2 - \ g» v E p % (3.5) 

and 

V Al V^g' 1/ + Q' A (V (A Q^ ) +/A,) = 0, (3.6) 

respectively. The relative size of the errors E\ v and f\ v is not clear. More 
important is that in this form the subsidiary system shows its non-linearity. 
If Q 1 ^ and dt Q' v vanish on S, equation l|3.6() will still imply that vanishes in 
the domain of dependence. If Q'^ and dt Q' v do not vanish on S, however, the 
following observations indicate that (|3.6(l may imply a growth of the solution 
which is worth than exponential. We assume g = dt 2 — 5 a \, dx a dx b (so that (|3.6|l 
decouples from l|3.3[0 . f\ v — 0, and initial data which satisfy for a = 1, 2, 3, 

Q' a =0, d t Q' a = 0, d a Q' a = 0, d a d t Q' o = on {t = 0}, 

so that the error resides only in the constant functions a = Q' , b — dtQ' Q on 
{' <>!• 

Equation (|3.6|l then implies Q' a = and reduces in fact to dt Q' Q = \ (c — 
Q 2 ) with c = 2 6 + a 2 . The integration gives Q' = a = const, if b = 0, 

v° = ^ a XTlt if * 25 > - fl2 ' Q ° = * if a2 = ~ 2h > and 



c\ tan ^ _ v l c l ' ^_ arc t an _a ^ if 26 < —a 2 . The solutions thus re- 



<\c\_ 

main bounded for t > if b > or if a > and —a 2 < 26 < 0, while they 
develop po/es at some > if 6 < and a < or if 26 < —a 2 and a > 0. 
(If the corresponding initial data would be modified outside the intersection of 
the hypersurface {t = 0} with the backward light cone of the point (i*, x a ), the 
solution would still become singular at (t^,x a ).) 

We note that in general any Killing vector field K of a vacuum solution g^ 
satisfies equation (|2.4|l with — in the coordinates x^. In the present 
case this gives solutions which grow linearly. 

For us the following observation is important: 

With our assumption on g M „ it follows that each neighbourhood of the initial 
data Q a — and dtQ a = contains initial data for \2.$ for which the solution 
Qaix^) become unbounded at some finite x° = t* > 0. 
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4 Concluding remarks 



The only purpose of the discussion above was to indicate the behaviour of solu- 
tions to the subsidiary system. We expect to find a similar singular behaviour 
of the solutions to equation (|2.4|) if the latter is given with a general smooth 
metric g^ v (extending sufficiently far into the future). It would be interesting 
to know whether the singular data set, i.e. the subset of data which determine 
singular solutions of equation (|2.4II , is open or of lower dimension in the set of 
all data if the metric g is given. While the singular data set itself will depend on 
the metric g, such a characterization may be independent of the chosen metric 
and may help avoid entering the singular sector. 

In a consistent discussion of the growth of one would have to consider the 
main system and the subsidiary system 

R» v = V (/1 Q U ) , V p V Q v + Q x V (A Q v) = 0, 

as a coupled system. The main equation can be studied independently. The 
subsidiary system, implicit in the main system, depends on the metric defined 
by the latter. If tends to grow, the main system will react to it, whether 
for the better or the worth is not clear. Whether with changing g^ v the non- 
linearity of the subsidiary system can still imply a blow up of at a finite 
time needs to be analysed. 

For this purpose it might be interesting to study under simplifying assumptions 
such as spherical symmetry whether the model system 

considered as Einstein equations coupled to a source field given by a vector field 
q^, will develop a blow up for suitable data. The reduction of these equations 
is obtained by a slight modification of the one described above. The system 
can be simplified further by assuming to be a differential — Vp / of some 
function /. The second equation will then be implied if / satisfies 

V M V/ + V M /V' 1 / = const., 

which is a wave map equation in the case where the constant on the right 
hand side vanishes. In any case the results might lead to an identifcation of 
a mechanism responsible for the growth of constraint violations and to the 
development of methods to avoid them. 

Some authors add terms built from Q v to the main system, which appear to 
reduce the growth of the constraint violations in certain calculation [3] , [7j . This 
may be related to the different effects of the non-linearities which are obtained 
in the appropriately modified subsidiary systems. 

There are available now many different types of reductions. Depending on 
several choices, the subsidiary system may control the preservation of the gauge 
or the preservation of constraints or a mixture thereof. In spite of the different 
appearance of the resulting main and subsidiary systems we expect that similar 
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non-linearities as the one discussed above will occur in any subsidiary system, 
though, depending on the system, they may have different effects. 
In a numerical scheme for the second order wave equations the subsidiary sys- 
tem, which is of third order in the metric, can, of course, hardly be identified 
any longer as a kind of identity and the relations between the two systems is 
obscured. But if the non-linearity of the subsidiary system can have for non- 
vanishing initial data and <9 t Q M drastic effects in the continuum model, they 
are likely to be reflected in numerical calculations. 
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